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Abstract 

Accurate total reaction cross section models are important to achieving reliable predictions from spallation and transport codes. 
The latest version of the Cascade Exciton Model (CEM) as incorporated in the code CEM03.03, and the Monte Carlo N-Particle 
transport code (MCNP6), both developed at Los Alamos National Laboratory (LANL), each use such cross sections. Having 
accurate total reaction cross section models in the intermediate energy region (~50 MeV to ~5 GeV) is very important for different 
applications, including analysis of space environments, use in medical physics, and accelerator design, to name just a few. The 
current inverse cross sections used in the preequilibrium and evaporation stages of CEM are based on the Dostrovsky et al. model, 
published in 1959. Better cross section models are available now. Implementing better cross section models in CEM and MCNP6 
should yield improved predictions for particle spectra and total production cross sections, among other results. Our current results 
indicate this is, in fact, the case. 
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1. Introduction 

Total reaction cross section models have a significant im¬ 
pact on the predictions and accuracy of spallation and trans¬ 
port codes. The latest version of the Cascade Exciton Model 
(CEM) mm as incorporated in the code CEM03.03 [|2][3], and 
the Monte Carlo N-Particle transport code (MCNP6) (4), both 
developed at Los Alamos National Laboratory (LANL), each 
use such cross sections for different purposes. While total reac¬ 
tion cross sections are used throughout the transport and spal¬ 
lation models, there are two main utilizations. MCNP6 uses 
total reaction cross sections to determine where a reaction oc¬ 
curs (through the mean-free path length), and then with what 
nucleus the projectile interacts with, and lastly what type of 
interaction it is (inelastic or elastic). CEM uses total reaction 
cross sections as inverse cross sections to predict what the ex¬ 
cited nucleus emits. Phenomenological approximations of total 
reaction cross sections are also used by CEM03.03 as the de¬ 
fault option for normalization of all results in the case of reac¬ 
tions induced by protons and neutrons, when CEM03.03 is used 
as a stand alone code, outside any transport codes; see details 
in Refs. (2l|3]|. 

Having accurate total reaction cross section models in the 
intermediate energy region (~50 MeV to ~5 GeV) is impor¬ 
tant for many different applications. Applications in space in¬ 
clude astronaut radiation dosage, electronics malfunction anal¬ 
ysis, structural materials analysis, and Galactic Cosmic Rays 
(GCRs) shielding. Medical applications include hadron ther¬ 
apy for cancer a , radiation shielding, medical isotope produc¬ 
tion, and high-radiation environment dosimetry. Other appli¬ 
cations include accelerator design and simulation. In addition. 
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implementing better inverse cross sections in CEM should pro¬ 
vide more reliable predictions; that is, our current work should 
be useful also from an academic point of view, allowing us to 
better understand the mechanisms of nuclear reactions. Lastly, 
the 2008-2010 IAEA Benchmark of Spallation Models recom¬ 
mended an improvement to CEMs ability to predict the pro¬ 
duction of energetic light fragments ®ia. Our improvement 
of the inverse cross sections used by CEM03.03 addresses di¬ 
rectly this point, both for a better description of light fragments, 
but also of nucleons. 

The current inverse cross sections used in the preequilibrium 
and evaporation stages of CEM are based on the Dostrovsky 
et al. model, published in 1959 [8). (For more information 
about the stages of CEM in its model of spallation reactions, 
see Ref. mmm .) Better total reaction ( inverse ) cross section 
models are available now COM]- 

MCNP6 uses an update of the Barashenkov and Polanski 
(B&P) cross section model lf2Tj as described briefly in Il22ll23l 
to calculate the mean-free path length for neutrons, protons, and 
light fragments up to 4 He. It uses a parameterization based on 
a geometric cross section for light fragments above 4 He. Im¬ 
plementing better cross section models in CEM and MCNP6 
should yield improved results of particle spectra and total pro¬ 
duction cross sections, among other results. Our current results, 
upgrading the inverse cross section model in the preequilibrium 
stage of CEM, prove that this is, in fact, the case. 

This cross section development work is part of a larger 
project aimed at enabling CEM to produce high-energy light 
fragments (9] [24] 25). Figs. |T| and [2] illustrate two examples of 
results of that project: comparing results from CEM03.03 with 
an upgraded Modified Exciton Model (MEM) to results from 
CEM03.03 unmodified. For some reactions we obtained good 
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results (see, e.g.. Fig. [TJ, and for other reactions, while our re¬ 
sults showed improvement, they could still be better (see, e.g.. 
Fig# We decided to upgrade the inverse cross section models 
used by CEM, in the preequilibrium stage, to improve such re¬ 
sults further. As CEM is the default event generator in MCNP6 
in the intermediate energy range, once these results our imple¬ 
mented into MCNP6 (to be completed soon), we should see a 
corresponding improvement in MCNP6 as well. 



Figure 1: Comparison of experimental data by Machner et al. on (green cir¬ 
cles) with results from the unmodified CEM03.03 (blue dotted lines) and the 
modified-MEM CEM03.03, we refer to here as CEM03.03F (24l | 25 ) (red solid 
lines) for 27 Al(p, 6 Li)X with incident proton energy of 200 MeV and for emis¬ 
sion angles of 20°, 45°, 60°, 90°, and 110°. 



Figure 2: Comparison of experimental data by Budzanowski et al. CD (green 
circles) with results from the unmodified CEM03.03 (blue dashed lines) and 
the modified-MEM CEM03.03, we refer to here as CEM03.03F f^4ll25l (red 
solid lines) for 197 Au(p, 7 Li)X with incident proton energy of 1200 MeV and for 
emission angles of 15.6°, 20°, 35°, 50°, 65°, 80°, and 100°. The 100°spectrum 
(the lower set) is shown unsealed, while the 80°, 65°, etc., down to 15.6°spectra 
are scaled up by successive factors of 10, respectively. 


2. Background 

As mentioned above, the current inverse cross sections in 
CEM are based on the Dostrovsky et al. model ED- It is based 
on the strong absorption model and its general form is as shown 
in Eq. #. 


o-Dost. = nrlA 2,i aj(l - y). ( 1 ) 

The Dostrovsky et al. model was not intended for use above 
about 50 MeV/nucleon, and is not very suitable for emission 
of fragments heavier than 4 He. Better total reaction cross sec¬ 
tion models are available today, most notably the NASA model 
mm. The NASA (or Tripathi et al.) model is also based 
on the strong absorption model and its general form is shown 
in Eq. 0- 5j, X m , and Bj are discussed more fully later, and 
defined in Eq. Q. The NASA cross section attempts to sim¬ 
ulate several quantum-mechanical effects, such as the optical 
potential for neutrons (with X m ) and collective effects like Pauli 
blocking (through 6t). (For more details, see Refs. | fT0UT2l .) 

o~nasa = + A X j 2 + djYd - R c ~— )X m , (2) 

' cm 

where 

/o is a constant related to the radius of a nucleus; 

Ap is the mass number of the projectile nucleus; 

Ap is the mass number of the target nucleus; 

6t is an energy-dependent parameter; 

R c is a system-dependent Coulomb multiplier; 

B T is the energy-dependent Coulomb barrier; 

T cm is the colliding system center-of-momentum energy; 

X,„ is an optical model multiplier used for neutron-induced 
reactions. 

There are other proposed total reaction cross section models, 
such as those by Shen, et al. ED- and Takechi, et al. El, 
amongst others HMD- It should be noted that both the Shen 
model and the Kox model have projectile-target assymetry, as 
discussed in Ref. (28) , In Ref. Il29l . Sihver et al. explores a new 
total reaction cross section used in PHITS: the hybrid Kurotama 
model. This model is a combination of the Black Sphere model 
071 and the NASA model IflOHTZI . Ref. |[30ll compares a num¬ 
ber of different total reaction cross section models, most notably 
those in FLUKA, NASA, and several other recently developed 
models. 

PHITS uses the NASA model as its default total reaction 
cross section model, but Shen can be specified as an option 
EOl. FLUKA uses a modified version of the NASA model as 
its total reaction cross section model ED- GEANT4 has the 
option to use NASA, or a number of other total reaction cross 
section models such as Shen El or Sihver l32l . or the Axen- 
Wellisch E3l total reaction cross section parameterizations for 
high-energy hadronic interactions. See Ref. E41 [35l for more 
details on the total reaction cross section models used in PHITS, 
FLUKA and GEANT4. 

Stepan Mashnik with collaborators l36l [37l and Dick Prael 
with coauthors l23l [381 previously conducted at LANL an ex¬ 
tensive comparison of the NASA IToUTZl . Tsang et al. fl9l . 
Dostrovsky et al. ED, Barashenkov and Polanski (using their 
code called CROSEC) IfZTl . and Kalbach l20l systematics for 
total reaction (inverse) cross sections. Fig. [3] illustrates some 
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Figure 3: Absorption (inverse) cross section by energy for various reactions, as calculated in Ref. 1361 by the NASA I10U121 , Kalbach (23, and Dostrovsky et al. 
Q3 systematics, as well as with a “Hybrid approach” suggested in [36) to account for both NASA 1111 and Kalbach 123 systematics, in the case of neutron-induced 
reactions. “BAR93” shows experimental data from Ref. OH; “DUB89” shows data from Ref. (40); and “AUC94” shows data from Ref. ED. 
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results from the study (36). Their results found that the NASA 
total reaction cross section model was superior, in general, to 
the other available models. See Ref. ED EH ED [42) for details 
of their findings. 

3. Comparison of Total Reaction Cross Section Models 

We built in CEM03.03F the NASA model Q0H3II and the 
models used in the preequilibrium (labeled as “Dostrovsky” in 
our figures below) and the evaporation (described with the code 
GEM2 by Furihata (431 , therefore labeled in our figures below 
as “GEM2”) stages of CEM03.03, and also compared reactions 
to calculations from the Barashenkov and Polanski (B&P) sys- 
tematics EL and, for comparison, to two neutron- and proton- 
induced reaction cross sections calculations by MCNP6 |4) . 
Note that MCNP6 uses currently an updated and improved ver¬ 
sion of the initial Barashenkov and Polanski (B&P) systematics 
EE as outlined briefly in Refs. ED EH, to simulate the mean- 
free path length of nucleons in matter. 

3.1. Neutron-Induced Reactions 

Fig. [^displays the total reaction cross section for n + 208 Pb, 
as calculated by the NASA, Dostrovsky et al., GEM2, and 
B&P models, and compared to calculations by MCNP6 and 
experimental data. There are several things to notice: 1) the 
Dostrovsky and GEM2 (also a Dostrovsky-based model) both 
approach asymptotic values very quickly-thus they are not 
as useful at their constant values, and 2) the NASA model, 
while much better at predicting the total reaction cross section 
throughout the energy region of projectiles, falls to zero at low 
energies in the case of neutrons, where we do not have Coulomb 
barriers. For this reason, we can not use the NASA model as 
an approximation for inverse cross sections in the case of low- 
energy neutrons: neutrons are emitted with low energies, too. 
In the case of low energy neutrons, we use the Kalbach system¬ 
atics [20], which proved to be a very good approximation for 
the inverse cross section of low-energy neutrons, as discussed 
in Ref. [36] and in sub-section 4.1 below. Note that this prob¬ 
lem of neutron cross sections was addressed first for the code 
CEM2k in Ref. (36) . by combining the NASA systematics by 
Tripathi, Cucinota, and Wilson EDHia and the Kalbach param¬ 
eterization (20) into a FORTRAN routine called hybrid. We ad¬ 
dress this problem here, for our current CEM03.03F code, in a 
very similar way (see Ref. (42) '). 

See Ref. (42) for results of other neutron-induced reactions. 

3.2. Proton-Induced Reactions 

Fig. 0 illustrates calculated total reaction cross sections by 
the NASA, Dostrovsky et al., GEM2, and B&P models, com¬ 
pared to calculations by MCNP6 and experimental data. The 
NASA model appears to be superior to the Dostrovsky-based 
models. 

As we can see from Figs. 5 and 4 on nucleons, as well 
as from examples on complex-particles and fragments heavier 
than 4 He presented below in Figs. 6 and 7, and in numerous 
figures published in Refs. [23, 36-38], the Barashenkov and 


n + 208 Pb 



Figure 4: Reaction cross section for n + 208 Pb, as calculated by the NASA, 
Dostrovsky et al., GEM2, and B&P models. The black dots are cross section 
calculations of MCNP6, and the yellow circles are experimental data Ena. 

Polanski approximations also agree very well with available 
data. For this reason, the B&P parametrization was chosen to 
be used for the calculation of the total reaction cross sections 
in the transport code MCNP6 [4], and in several other transport 
codes, too, as far as we know. Howerver, our numerous current 
comparisions for various reactions, as well as the voluminous 
results published in Refs. [23, 36-38], show that, on the whole, 
the NASA approximation agree a litlle better with most of the 
available experimental data than the B&P systematics does. 

See Ref. (42) for results of other proton-induced reactions. 

3.3. Heavy-Ion Induced Reactions 

We never tested before how CEM03.03 calculates inverse 
cross sections for light fragments (LF) heavier than 4 He. We 
address this question below. 

Fig. 0 illustrates calculated total reaction cross sections by 
the NASA, Dostrovsky et al., GEM2, and B&P models for the 
reactions a + 28 Si and 6 Li + 208 Pb, compared to experimental 
data. 

Fig. [^displays the total reaction cross section for 12 C + l2 C, 
as calculated by the NASA, Dostrovsky et al., GEM2, and B&P 
models and compared to experimental data and to measured to¬ 
tal charge-changing (TCC) cross sections. TCC cross sections 
should be 5% - 10% less than total reaction cross sections, as 
TCC cross sections do not include the neutron removal cross 
section. 

See Ref. (42) for results of other heavy-ion-induced reac¬ 
tions. 
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Figure 5: Reaction cross section for p + 12 C, as calculated by the NASA, 
Dostrovsky et al, GEM2, and B&P models. The black dots are cross section 
calculations of MCNP6, and the yellow circles are experimental data ED 

We determined that the NASA cross section model fits the 
experimentally measured data, in general, better than the other 
models tested. 


4. Implementation of NASA Cross Section Model into 

CEM03.03F 

The implementation of the NASA cross section model into 
CEM involved adding Kalbach systematics for low-energy neu¬ 
trons, updating the emission width calculation, and upgrading 
the emitted fragment kinetic energy simulation. 

4.1. Kalbach Systematics 

We added in CEM03.03F the Kalbach systematics |[20l to re¬ 
place the NASA inverse cross sections mm for low-energy 
neutrons, similar to what was suggested and done in Ref. | [36] 
for the code CEM2k. Fig. [8] displays the Kalbach systemat¬ 
ics implementation for the cross section n + 208 Pb. At around 
24 MeV and below, the calculation switches to Kalbach sys¬ 
tematics, and uses the NASA model throughout the rest of the 
neutron-energy range. The Kalbach systematics is scaled to 
match the NASA model results at the switchpoint so as not to 
have a large jump. 

As part of the Kalbach systematics implementation in 
CEM03.03F, switchpoints and scaling factors must be obtained 
for all possible residual nuclei, by mass number. Ref. Il42ll pro¬ 
vides tables of these. 


a + 28 Si 



6 Li + 208 Pb 



Figure 6: Reaction cross section for a + 28 Si and 6 Li + 208 Pb, as calculated by 
the NASA, Dostrovsky et al., GEM2, and B&P models. The yellow circles are 
experimental data 1564601 . 

4.2. Emission Width, T j, Calculation 

CEM uses the inverse cross section, crj 1 ’, in determining 
what particles and/or fragments are emitted from the excited 
nucleus. We use the total reaction cross section as the best ap¬ 
proximation for this inverse cross section. The emission width 
Fj, or the probability of emitting fragment type j, is calculated 
according to Eq. d3|. It is dependent upon cr"' v (see more details 
in Refs. U 
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Figure 7: Reaction cross section for 12 C + 12 C, as calculated by the NASA, 
Dostrovsky et al., GEM2, and B&P models. The yellow circles are experimen¬ 
tal data Q3H23HD and the blue squares are total charge-changing cross section 
(TCC) measurements 16211631 . 


n + 208 Pb 



Figure 8: Reaction cross section for n + 208 Pb, as calculated by the NASA, 
NASA-Kalbach hybrid (black line), Dostrovsky et al ., GEM2, and B&P mod¬ 
els, as indicated. The black dots are cross section calculations by MCNP6, and 
the yellow circles are experimental data from Refs. 


r - H i 2sj + i M 

(p, h,E) = J pj^(p, h) 

oXp- l,h,E-Bj-T ) 
u)(p , h, E) 


(3) 


To-‘" v (T)dT , 


where 


sj is the spin of the emitted particle j; 

/jj is the reduced mass of the emitted particle /; 
oj is the level density of the n-exciton state; 

Bj is the binding energy; 

Vj is the Coulomb barrier; 

T is the kinetic energy of the emitted particle j; 
crj v is the inverse cross section; and 

% creates zero probability of emission if the number of parti¬ 
cle excitons is less than the number of nucleons of parti¬ 
cle j. 

Eq. ([3ji is written in its simplest form, as is valid for neutrons 
and protons only. An extension of Eq. 0 for the case of com¬ 
plex particles and light fragments (LF) is described in detail in 
Ref. Q. 

In the “standard” (i.e., “old,” for this study) calculation by 
CEM03.03, performed with a FORTRAN routine called gam- 
agu2, therefore referred to below as “gamagu2,” the Dostro¬ 
vsky et al. form of the inverse cross section is simple enough 
so that for neutrons and protons this integral can be done an¬ 
alytically. However, for complex particles, the level density, 
to , becomes too complicated (see details in Refs. |QJ|3]), there¬ 
fore the integral is evaluated numerically. In this case, a 6-point 


Gaussian quadrature is used when the exciton number is 15 or 
less, and a 6-point Gauss-Laguerre quadrature is used when the 
number of excitons is over 15. We will soon see why the two 
different integration methods are needed. 

In our current calculations we adopt here for CEM03.03F 
(performed with a FORTRAN routine called gamagu3, here¬ 
after referred to as “gamagu3”), the NASA form of the cross 
section is too complicated and the integral is always calculated 
numerically. We use an 8-point Gaussian quadrature when the 
number of excitons is 15 or less, and an 8-point Gauss-Faguerre 
quadrature when the number of excitons is greater than 15. 

The partial transmission probability Aj, or the probability that 
a particle or a fragment of the type j will be emitted with kinetic 
energy T, is equal to the integrand of Eq. #• For the emission 
of FF this is equal to 


Aj(p,h,E,T ) -y } pj%(p,h) 

to(p - pj, h, E - Bj - T) 


u(Pj 


to(p, h, E) 
, 0, T + Bj) 


Sj 


To-'P’iT ), 


(4) 


where 


V(2/i.-) 3/2 ... 

^=-^~ (2s J + l)(T + B J ) ■ (5) 

See Ref. Il64l for details on Eq. d5l. For completeness sake, 
we write Eq. (|4]) (and also Eq. (j6jT below in the text) in their 
“complete form,” as they should look in the case of complex 
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particles and LF, but not in their “simplest” version needed only 
for nucleons as exemplified by Eq. <0 

As an example, Fig.[9]shows Aj for the emission of neutrons 
from a 198 Au excited nucleus, with an internal nucleus energy 
U of 200 MeV, using either the Dostrovsky et al. or NASA 
cross section. The top plot is for 55 excitons and the bottom 
plot is for 10 excitons. Notice that for high exciton number, A j 
becomes more concentrated in the low-energy region. Table [T] 
displays the abscissas for an 8-point Gaussian and an 8-point 
Gauss-Laguerre quadrature. For a small number of excitons 
(< 15) the Gaussian quadrature performs adequately. However, 
we see that in the 55-exciton case the Aj becomes negligible 
by about 30 MeV, requiring a different integration method. For 
high-exciton number the Gauss-Laguerre integration method is 
a much better choice than the simple Gaussian. 


Table 1: 8-point Gaussian and Gauss-Laguerre sampling points 


8-pt Gaussian 

8-pt Gauss-Laguerre 

3.84 MeV 

0.428 MeV 

19.7 MeV 

2.27 MeV 

45.9 MeV 

5.66 MeV 

79.0 MeV 

10.7 MeV 

114. MeV 

17.7 MeV 

148. MeV 

27.1 MeV 

174. MeV 

39.6 MeV 

190. MeV 

57.5 MeV 


Fig. 10 shows a comparison of the simple Gaussian and 
Gauss-Laguerre quadratures for 55 excitons. This figure also 
displays Aj for the NASA-Kalbach cross section. Notice that 
the NASA-Kalbach has much higher values of A„ at the low 
end of the spectrum than the pure NASA. The purple dots are 
the 8-pt Gaussian quadrature and the black dots are the 8-pt 
Gauss-Laguerre quadrature. The Gaussian was exceptionally 
fortunate in that it struck the peak with its one low-end point. 
However, this leads to significant overestimation of Aj down 
the tail. The Gauss-Laguerre underestimates the peak but then 
overestimates slightly along the tail. Even though it is clear 
this is not a very close fitting of Aj, changing to a 10-pt Gauss- 
Laguerre only yielded a 0.2% difference. A future project could 
include investigating the behavior of Aj across the variable land¬ 
scape, and implementing an adaptive quadrature scheme. How¬ 
ever, whatever numerical integration method we use, it must be 
fast as this integral is calculated hundreds of times for every 
event, and therefore billions of times for a typical simulation. 

Fig. 11 shows the plots of T,- as a function of the internal en¬ 
ergy of the excited nucleus for emitted neutrons, protons, and 
4 He from an excited 198 Au nucleus with 55 excitons, 25 particle 
excitons, and 13 charged particle excitons. “Gamagu2” shows 
the old CEM03.03 Vj calculation results. “Gamagu3” shows 
the results of our new calculations, using either the Dostrovsky 
et al. or NASA inverse cross sections. Note that “Gamagu2” 
should be very similar to “Gamagu3-Dostrovsky” because the 
only significant difference is the method of integration. The 
proton and neutron Tj differences between “Gamagu2” and 
“Gamagu3-Dostrovsky” arise from numerical integration used 




Figure 9: Aj as a function of the kinetic energy of the emitted neutron, from an 
excited 158 Au nucleus with U = 200 MeV and 55 excitons, 25 particle excitons, 
and 13 charged particle excitons (top plot) and 10 excitons, 6 particle excitons, 
and 3 charged particle excitons (bottom plot). 


in our new FORTRAN routine gamagu3 of CEM03.03F ver¬ 
sus an analytical calculation used in the CEM03.03 FORTRAN 
routine gamagu2. 

Better integration methods could be investigated at a later 
time. However, current integration methods are sufficient be¬ 
cause individual I j precision is not extremely important for 
choosing what type of particle/LF j will be emitted. In contrast 
to analytical preequilibrium models, the Monte Carlo method 
employed by our CEM uses the ratios of F 7 - to the sum of F ; 
over all j. That is, if we estimate all T ; with the same percentage 
error, the final choice of the type j of particle/LF to be emitted 
as simulated by CEM would be the same as if we would cal- 
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Figure 10: Aj as a function of the kinetic energy of the emitted neutron, from an 
excited 198 Au nucleus with U = 200 MeV and 55 excitons, 25 particle excitons, 
and 13 charged particle excitons. 


culate all exactly. We think that this is the main reason why 
CEM provided quite reasonable results using the old Dostro- 
vsky et al. approximation for inverse cross sections, in spite of 
the fact that, as we see from Figs. m individual inverse cross 
sections calculated with the Dostrovsky et al. method are not 
good enough in a large range of energies. The ratios Fj/ E/T/) 
were probably estimated well enough, providing a reasonable 
Monte Carlo sampling of j. 


4.3. Kinetic Energy Simulation 


Once a fragment type j has been chosen for emission, the 
kinetic energy of this fragment needs to be determined. This is 
done by sampling the kinetic energy from the Aj distribution. 
Our new Aj, with the NASA cross section, is: 


Aj(p, h, E, T) =yj S J 2 ^ 3 p,%(p, h) 

cj(p - pj , h, E - Bj - T) 


co(p, h, E) 
ca(pj,0,T + Bj) 


( 6 ) 


Titri 


gj 


x (A 1 ' 3 + A 1 ' 3 + 6 t ) z ( 1 -R c -Y-)X m {T). 


11/3 


B T 





U (MeV) 

Figure 11: as a function of the internal energy of the excited nucleus for 

emitted neutrons, protons, and 4 He from an excited 198 Au nucleus with 55 ex¬ 
citons, 25 particle excitons, and 13 charged particle excitons. 


where gj is defined by Eq. (|5j» and 
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6 t = 1.855 + - D[ 1 - e~ TIT '] - 0.292e~ Tp92 


-1/3 


0.91 (A 7 - - 2Z T )Z P 
A T A P 


x m = 1 - Xi exp 


(7) 


B t = \MZpZ T l 


Xi(1.2+ 1.6[1 -exp(T/15)]))’ 
l.2(A l J 3 + A\! 3 )) 


Yp + Yp + 


-1/3 


The details of ro,R c ,S,D, T\,rp, and r T can be found in lfT2l . 
Note that the NASA inverse cross sections contain depen¬ 
dences on both the lab-reference-frame kinetic energy (T) and 
the center-of-momentum-reference-frame kinetic energy (T, 
The relativistic transformation between the two is not trivial. 
In addition, T is in units of MeV/nucleon in the NASA in¬ 
verse cross sections, while T cm is in units of MeV. The level 
density, to, also contains T -dependences, also in units of MeV. 
Finally, as noted above, for neutrons we use a NASA-Kalbach 
(“hybrid”) inverse cross section in place of the pure NASA ap¬ 
proximation. To conclude, the energy-dependence of Aj for 
our new NASA-Kalbach inverse cross section approximation is 
very complicated, which affects the method we chose to sample 
Tj, as discussed below. 

To sample Tj uniformly from the Aj distribution using the 
Monte Carlo method, we must first find the maximum of Aj. 
In CEM03.03, this is done analytically using the derivative of 
Aj with respect to T due to the simple nature of the energy- 
dependence in the systematics by Dostrovsky et al.. As pre¬ 
viously explained, however, the NASA cross section energy- 
dependence is extremely complicated and therefore we find the 
maximum of Aj numerically using the Golden Section method. 
This also provides us flexibility in the future to modify Aj with¬ 
out consequence to our kinetic energy module. 

After finding the maximum value of Aj, the kinetic energy of 
the emitted fragment j is uniformly sampled from the Aj distri¬ 
bution using a Gamma distribution (shape parameter a = 2) as 
the comparison function. Fig.[l2]illustrates results for the prob¬ 
ability of emitting 6 Li with a given kinetic energy 7) ,. Proba¬ 
bilities from the Aj distributions with the NASA inverse cross 
sections differ slightly from those with the Dostrovsky et al. 
inverse cross sections primarily because the NASA coulomb 
barriers are based on T cm , as opposed to T. 


5. Results 

Our preliminary results are promising. Fig. [13] shows the 
double differential cross section for the production of 6 Li and 
7 Be from the reaction 200 MeV p + 59 Co. Notice the improved 
agreement with data in the high-energy tails. This reaction also 
highlights the importance of eventually upgrading the inverse 
cross sections used in the evaporation stage of CEM as well. 
The evaporation stage produces the peak of the spectra, which 
for this reaction is too low, especially for 7 Be. With the im¬ 
plementation of the NASA inverse cross sections in the pree¬ 
quilibrium stage we see improved agreement with data in the 



Figure 12: The normalized probability of emitting 6 Li with a given kinetic en¬ 
ergy 7/,, simulated using the Monte Carlo method according to Eq.|6]in the 
preequilibrium stage. The circles are results from the old kinetic energy sub¬ 
routine; the squares are results from the new kinetic energy subroutine using the 
Dostrovsky et al. inverse cross section; the triangles are from the new kinetic 
energy subroutine using the NASA inverse cross section. 


high-energy tails, but in order to achieve improved agreement 
in the peak we would need to also implement the NASA inverse 
cross sections in the evaporation stage. We plan to do this in the 
future. 

For another example of our results, Fig.[l4]displays the dou¬ 
ble differential cross section for the production of 6 He and 7 Li 
from the reaction 1200 MeV p + 197 Au. The blue dashed 
lines are the expanded-MEM (i.e., CEM03.03F) results with the 
Dostrovsky et al. inverse cross sections, and the red solid lines 
are results by CEM03.03F with the upgraded NASA-Kalbach 
(i.e., “hybrid”) inverse cross sections. The green circles are ex¬ 
perimental data from Ref. |[27l . We see an improved accuracy 
in the high-energy tails of spectra calculated with the NASA in¬ 
verse cross sections, although some of the results are too hard 
and there is a dip in the spectra at 50-75 MeV. We would like 
to note a recent paper by A. Boudard et al. Il65ll . which obtained 
very similar results for 7 Li using INCL4.6 + ABLA07, and sim¬ 
ilar results for 6 He but with a little lower evaporation peak. Cur¬ 
rent work is being undertaken to expand the coalescence model 
in CEM, which helps soften the spectra and smooth out the dip 
that appears around 50-75 MeV. This work is an ongoing pro¬ 
cess, but we display some of our preliminary results in Fig. [15] 
For further details, see Refs. 11671 16811 . 
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200 MeV p + 59 Co -> ... 



Figure 13: Double differential cross section for the production of 6 Li and 7 Be from the reaction 200 MeV p + 59 Co, for the angles of 20°, 45°, 60°, 90°, and 110°. 
The 110° spectra (the lower sets) are shown unsealed, while the 90°, 60°, 45°, and 20° spectra are scaled up by successive factors of 10, respectively. The blue dashed 
lines are the expanded-MEM (i.e., CEM03.03F) results with the Dostrovsky et al. inverse cross sections, and the red solid lines are results by CEM03.03F with the 
upgraded NASA-Kalbach (i.e., “hybrid”) inverse cross sections. The green circles are experimental data by Machner, et al HD 


1200 MeV p + 197 Au -> ... 



Figure 14: Double differential cross section for the production of 6 He and 7 Li from the reaction 1200 MeV p + 197 Au, for the angles of 15.6°, 20°, 35°, 50°, 
65°, 80°, and 100°. The 100°spectra (the lower sets) are shown unsealed, while the 80°, 65°, etc., down to 15.6°spectra are scaled up by successive factors of 10, 
respectively. The blue dashed lines are the expanded-MEM (i.e., CEM03.03F) results with the Dostrovsky et al. inverse cross sections, and the red solid lines are 
results by CEM03.03F with the upgraded NASA-Kalbach (i.e., “hybrid”) inverse cross sections. The green circles are experimental data from Ref. EzJ- 


6. Conclusion 

The inverse cross section approximation in the preequilib¬ 
rium and evaporation stages of CEM03.03 is based on the 
Dostrovsky et al. inverse cross section model. Better cross sec¬ 
tion systematics are available at present. We performed a com¬ 


parison of several inverse cross section models and determined 
that the NASA (Tripathi, et al.) approximation is, in general, 
the most accurate when compared with experimental data. 

We implemented the NASA inverse cross section model into 
the preequilibrium stage of CEM03.03F. This included writing 
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sections. 



Figure 15: Comparison of experimental results of the reaction 480 MeV p + 
nar Ag —> 6 Li at 60° by Green et al. EH (green circles), with simulations from 
the original CEM03.03 (brown dashed-dotted lines), CEM03.03F without co¬ 
alescence expansion (blue solid lines) and the CEM03.03F with coalescence 
expansion (red dashed lines). 

FORTRAN modules containing the NASA total reaction cross 
section and coulomb barrier approximations, adding Kalbach 
systematics for low-energy neutron inverse cross sections, re¬ 
writing the r ; routines (including transforming them into mod¬ 
ular FORTRAN), adding Gauss-Laguerre quadrature for cases 
of high exciton number, and re-modeling the selection of parti¬ 
cle or light fragment kinetic energy. These technical improve¬ 
ments lead to greater flexibility and robustness, and future up¬ 
grades can be made easily. 

Our preliminary results are promising and indicate improved 
agreement with experimental data using the NASA inverse 
cross section model versus the Dostrovsky el al. approxima¬ 
tion. 

There are several implications of this work on MCNP6. 
CEM03.03 is the default event-generator in MCNP6 for high- 
energy collisions induced by nucleons, pions, and gammas at 
energies up to several GeVs. Improvements to the CEM inverse 
cross sections should, therefore, result in improved predictions 
of particle spectra and total production cross sections, espe¬ 
cially above ~100 MeV and for fragments heavier than 4 He, 
among other results. 

MCNP6 uses the updated Barashenkov and Polanski total re¬ 
action cross section systematics to simulate the mean-free path 
of neutrons, protons, and light fragments up to 4 He. It uses a pa¬ 
rameterization based on a geometric cross section for fragments 
heavier than 4 He. Possible direct improvement of MCNP6 may 
be obtained by replacing the Barashenkov and Polanski model 
with NASA systematics and by replacing the geometric cross 
section approach with the better NASA model. We hope to do 
this in the future. 

Future recommendations include investigating adaptive 
quadrature and upgrading the inverse cross section model used 
in the evaporation stage to the NASA-Kalbach (hybrid) cross 
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